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Abstract 



'■^ ■ Two alternative Vinen's evolution equations for the vortex line density L in counter- 

G I flow superfluid turbulence, are physically admissible and lead to analogous results in steady 

states. In Phys. Rev. B, 69, 094513 (2004) the most used of them was generalized to coun- 
terflow superfluid turbulence in rotating containers. Here, the analogous generalization for 
the alternative Vinen's equation is proposed. Both generalized Vinen's equations are com- 
pared with the experimental results, not only in steady-states but also in some unsteady 
situations. From this analysis follows that the solutions of the alternative Vinen's equation 

^SJ ! tend significantly faster to the corresponding final steady state values than the solutions 

t^^ ■ of the usual Vinen's equation, and that the latter seems more suitable for the description 

^^ I of the experimental available data. 



t^ ! 1 Introduction 

O 

There is experimental evidence that turbulent helium H has a peculiar behavior due to quan- 
K> i turn effects at some length scales, while at other length scales it appears similar to classical 

H ' hydrodynamic turbulence [I]-[5]. Quantum turbulence is described as a chaotic motion of 

quantized vortices in a disordered tangle. The measurements of vortex lines are described 
in terms of a macroscopic average of the vortex line length per unit volume L (briefly called 
vortex line density and which has dimensions length^^). 

The evolution equation for L under constant values of the counterflow velocity V (V =< 
Vn — Vg >, Vn and Vs being the velocities of the normal and superfluid components) and in 
absence of rotation was formulated by Vinen. Neglecting the influence of the walls, such an 
equation is [6]: 

^ = aVL^/^-/3KL^ (1.1) 

at 
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with V = |V|, K = h/m the quantum of vorticity {m the mass of the "^He atom and h Planck's 
constant, k ~ 9.97 10~^cm^/s) and a and /3 dimensionless parameters. 

Note however that another version of (jl.ip is the so-cahed alternative Vinen's equation, 
which is also admissible on dimensional grounds [7J-[9J: 

^ = A^^L-(3kL\ (1.2) 

at K 

The steady state solutions of (fTT]) and (fOI) are L = {oV/fiKf and L = AiV'^/^k^ 
respectively, in agreement with the experimental results in completely developed turbulent 
regime, which lead to L^'^ = jV/k, with 7 a dimensionless coefficient which depends on 
the temperature. Therefore, the difference between (jl.ip and (jl.2p must be searched in the 
dynamical aspects. This was carried out by Vinen himself (see sections 6 and 7 of Ref. [7]) and 
in more detail by Nemirovskii et al. [9] without arriving to definite conclusions, because the 
predictions of (II. ip and (II. 2|) in the domain of available experimental results are very similar 
to each other. 

Here we will look for a more general situation where the difference between p.ip and (|1.2p 
becomes enhanced. Essentially, equations (jl.ip and (jl.2p may be obtained from a microscopic 
approach based on vortex dynamics, though the derivation of the first one is more direct and 
straightforward than the other. Before proceeding, let us briefly recall that from a microscopic 
approach based on vortex dynamics the production term may be shown to be proportional 
to VL/Rpeak, where Rpeak is the intervortex spacing [51 [10]. Since the latter is of the form 
L^^/^, the form VL^'"^ adopted in (jl.ip follows in an immediate way. However, it has been 
also argued that [Sj, since in the steady state L^'^ is proportional to V, one could also write 
Rpeak 3-8 inversely proportional to V, in which case one would have for the production term 
the alternative form V^L adopted in (jl.2p . The equation (jl.2p was also derived by Lipniacki 
(pag. 177 of Refs. [2111]) through an alternative microscopic approach, which is referred to 
the reconnections of vortex lines. 

Let us say, in support of the form V'^L, that though it is less natural from a kinematical 
point of view, it is rather natural from an energetic point of view, if one assumes in agreement 
with the microscopical model^'^ that the friction force between the normal fluid and the vortex 
tangle is proportional to VL. In this case, the power delivered to the tangle per unit volume 
would be proportional to V^L. Since the energy density of the tangle is proportional to L, the 
production contribution to dL/dt would be proportional to the power delivered to the tangle, 
i.e. to V^L. 

Another motivation in support of (jl.2p is that it is closer to the phenomenological theory 
of classical turbulence [12] than equation (jl.ip . Indeed, by assuming that turbulence can be 
characterized by the line density L and that its derivative dL/dt is an analytic function of L, the 
relation AiV^L can be interpreted as the first term in a series expansion. However, equation 
(jl.ip has been much more used than equation (jl.2p . Both (jl.ip and (jl.2p are particular cases 
of: 

in fact for n = 1 one obtains (jl.ip . for n = 2 one has (1.2), and the cases 1 < n < 2 could 
correspond to fractal (intermittent) turbulence |13j-[15j. 

In recent years there has been growing attention in superfluid turbulence in rotating con- 
tainers 116] -[18], in which the formation of vortex lines is due both to the counterflow and 



the rotation, which has fostered the extension of Vinen's ideas to a wider range of situations 
1151119^ 120], In Ref. [19], a phenomenological generahzation of Vinen's equation (1.1) has been 
proposed for the evolution of L in the simultaneous presence of V and fi (fi being the angular 
velocity of the container) . A thermodynamic analysis to determine possible coupling terms 
between the evolution equations of L and V has been performed in Refs. [15]. 

Here, the extension of the form (II. 2p of Vinen's equation to rotating counterflow turbulence 
will be studied in order to explore whether this more general situation may provide further 
arguments to decide which of both starting equations, (II. ip or ()1.2p . is more suitable to 
describe actual experimental results. Of course, the final version will be fully satisfactory 
only when the macroscopic expression will be understood from a microscopical basis, and the 
respective coefficients of all terms will be microscopically calculated and found to coincide with 
macroscopic observations. However, this situation is still far ahead of our present abilities, 
because of the difficulties in modeling — in a statistically significative way — a system of 
vortices under rotation. Thus, a combined effort in macroscopic and microscopic perspectives 
seems a reasonable and promising way to proceed. 

The plan of the paper is the following. In Section 2 a new equation for the evolution of L in 
counterflow in rotating containers is written, through a modiflcation of the Vinen's alternative 
equation, and it is solved in steady and unsteady situations in order to compare it with the 
generalization of the usual Vinen's one made in Ref. [19]. In Section 3 a thermodynamic anal- 
ysis of counterflow rotational superfluid turbulence is performed, according to the formalism 
of nonequilibrium thermodynamics, to obtain the general form of the friction exerted by the 
tangle on the motion of the fluid. 

2 New equation for the dynamics of L in rotating counterflow 
superfluid turbulence 

There are not many experiments on counterflow in rotating containers. In the work of Swanson 
et al. [22], the counterflow velocity V was parallel to the rotation axis and the experimental 
observations consisted in measuring the attenuation of second sound, when it is propagated 
orthogonal to the rotation axis. They interpreted their results as measurements of the vortex 
line density L, and compared the observed line density with what would be expected if the 
two sources of vorticity (rotation and counterflow) simply added. Their results showed an 
interesting interplay between the ordered vortices of rotation and the disordered ones of coun- 
terflow. They found a critical counterflow velocity which marks the transition to a disordered 
turbulent tangle. In the limit of high rotation this critical velocity appears proportional to 

Another experiment is that carried out by Yarmchuck and Glaberson ^23j, in which a pair 
of horizontal parallel glass plates are arranged to form a large aspect ratio channel closed at one 
end with a heater nearby, and open at the other end to the liquid helium bath. The channel is 
rotated about a vertical axis orthogonal to the heat flux. In this way the counterflow velocity 
V is orthogonal to angular velocity il. By determining gradient of temperature and chemical 
potential as functions of V and J7, they found a linear regime, in which these gradients grow 
proportionally to the rotation speed, and a critical counterflow velocity associated with the 
onset of turbulent regime, which increases as V" when gets large. 



2.1 The evolution equation 

The mentioned experiments by Swanson, Barenghi and Donnelly [52] and by Yarmchuck and 
Glaberson [23] show that, when the heat flux and the rotation are simultaneously present, 
there appears a complex interaction between both processes in the formation and destruction 
of vortices. In the experiment performed in Ref. [22] (fi and V collinear), they observed 
that the effects of V and fi are not additive: in fact, for low values of O, the laminar regime 
(vortex- free regime) is absent and the total vortex line density is higher than L/j + Lh , Lr and 
Lh being the values of L in steady rotation and in steady counterflow superfluid turbulence 
respectively: 

20 9^2 

with 7 a dimensionless coefficient while, for higher values of Vt and V , the measured values of 
L are always less than Lh + L/j, and the deviation increases with V and $7. Therefore, the 
rotation facilitates the vortex formation, in the absence or for small counterflow velocities, but 
it hinders their lengthening for high values of V and Vt. 

For fast enough rotation, there are two critical counterflow-rotation velocities Vd and Vc2-, 
which scale as ^^/^ {Vd = Ci^/O, Vc2 = €2^^, with Ci = 0.053 cm sec'^/^, C2 = 0.118 
cm sec~^'^). For V < Vd, the length L per unit volume of the vortex lines is independent of 

V and agrees with the first expression in (|2.ip . For Vd ^V< Vc2, L is still independent of V 
and proportional to Q, with a slightly different proportionality constant than in the previous 
situation; finally, for V > Vc2, L increases and becomes proportional to V'^ at high values of 
V. 

Swanson et al. [22] interpreted the first transition as the Donnelly-Glaberson instability 
CCD,OG: excitation of helical waves (Kelvin waves) by the counterflow on the vortex lines 
induced by rotation, and the second as a transition to a turbulent disordered tangle. Tsubota 
et al. [18j also have paid attention to this experiment. They proposed that the regime Vd < 

V < Vc2 is a state of polarized turbulence, while for V > Vc2 the polarization is decreased by 
the large number of reconnections. 

In the regime of high rotation (0.2 Hz < J7/27r < 1.0 Hz and < F^ < 0.2 cm^/s^) and Jl 
parallel to V, equation (jl.ip has been generalized to incorporate the presence of rotation, as 
[191: 



^^ n t2 

—— = —pKL + ai 
at 



rl/2 V^ 

L ' — mi^-= 

K 



VL + P2 



rl/2 VO 



^L, (2.2) 



where mi and m2 are linked to the coefficients introduced in Ref. [19j by the relations mi = 
Pi/ai and m.2 = /3i//32, with the coefficients /?, ai, /3i, P2 and Pi depending on the polarization 
of the tangle, which was supposed function of and V. 

Note that, as production terms in Eq. (j2.2p a term in V and a term in i7^'^ were used; 
this was motivated by the dependence of the steady-state values of L^'^, in counterflow only 
and in rotation only, on V and on Jl^'^ (see equations ()2.ip ). and by the observation that the 
microscopic mechanism responsible for the growth of vortices (the mutual friction force) is the 
same in rotating helium II and in superfluid turbulence. There are present three destruction 
contributions: a term —(3kL'^ independent on V and on i7, present also in (11. ip and (II. 2p (this 
term, responsible for the vortex decay in pure counterflow, was determined by Vinen in analogy 
with classical turbulence) and the two terms — /JiilL and —(3^ ^^- L, which take account of 



the interactions between counterflow and rotation, reducing the length of the vortices; a term 
quadratic in V was neglected, because the values of V used in the experiments by Swanson et 
al. were not very high. 

The Eq. (j2.2p describes, in good agreement with experimental results, some of the most 
relevant effects observed in the experiments of Ref. [22]. However, as we have mentioned in 
the Introduction, the alternative Vinen's equation also describes well the experimental results 
in pure counterflow and therefore it is natural to ask how does it work when extended to 
incorporate rotation. 

Here, we suggest a new evolution equation for the evolution of vortex line density L in 
rotating counterflow, starting from the alternative Vinen's equation and following the lines of 
thought outlined in Ref. [19]. We consider the case in which V and are parallel to each 
other. The proposed new equation, reducing to (1.2) for vanishing rotation, is: 



at 
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n, (2.3) 



where the coefficients /3, Ai, z^i, Bi and 1^2 depend on the polarization of the tangle, which 
was supposed function of Q and V. 

We outline now a possible physical interpretation for the terms of production and de- 
struction of vortices introduced in this equation. The two production terms, Ai^L and 
BiQL, indicate that both rotation and counterflow favor the vortex formation: the quantities 
tn = k/{AiV'^) and tR = l/{BiQ) can be interpreted as the characteristic times for the for- 
mation of vortex lines due to the counterflow and to the rotation, respectively. As in (|2.2|) . 
three destruction terms are present. The term A\V\^V'^ j k} describes the complex interac- 
tion between rotation — which tends to straighten out the otherwise irregular vortex lines of 
the tangle, thus shortening them and reducing L — and counterflow, which randomize them. 
Another aspect especially worth of comment is the meaning of the two destruction terms, 
independent on y, —B\V2^li^ and —^kI?. One could argue, indeed, that at steady pure 
rotation there is no vortex destruction. Thus, in purely rotation situations the vortices are 
usually produced on the walls and they migrate to the bulk of the fluid in the cylinder; in this 
case, these terms would represent a repulsion force between parallel vortices, putting an upper 
limit to the possible number of straight vortices in the vortex array. 

In a general situation, these two destruction terms will incorporate real destruction of vor- 
tices due to breaking recombination of nonparallel vortices, and to repulsion between parallel 
segments of vortices in the presence of rotation. The superposition of these two different ef- 
fects is one of the reasons that the coefficients in the terms in l? and in O^ depend on the 
polarization of the tangle. The destruction term —[3kL'^, which appears also in (II. Ih . ()1.2p 
and ()2.2p . is not modified, in agreement with recent studies which show that the decay of the 
turbulence, in the absence of rotation and counterflow velocity, is analogous to that of classical 
turbulence [12] . 

As it was shown in Ref. [H], the simplicity of ()2.2p as (12. 3p is a little bit deceptive, because 
the coefficients appearing in it depend on the polarity and the anisotropy of the vortex tangle, 
which are taken as independent variables in the more detailed approach proposed in the recent 
paper [13]. Of course, this difficulty arises not only in the macroscopic approach, but also in 
the microscopic approaches. Since the coefficients appearing in (j2.3p depend on the anisotropy 
and on the polarization of the tangle, they have to depend on the angular velocity fi and on 
the counterflow velocity V. In particular, when = Eq. (j2.3p reduces to (1.2), and the 



coefficient Ai/(3 assumes the value: 



Ai 

13 



0.156. 



(2.4) 



In tfie analysis carried out in tfie present paper, the dependence of the coefficients on the 
polarity plays not an important role, because in Section 2.2 we are comparing two equations 
in a regime of values of V and Q (0.2 Hz < 0/27r < 1.0 Hz and < V'^ < 0.2 cm^/s^) with 
approximately the same polarity, and in Section 2.4 we compare the non-stationary behavior 
of the perturbations to a given physical situation, as described by two different equations. 

In the successive subsections we solve equation (12. 3|) in steady and unsteady situations and 
we show that it allows us to account for the experiments described in Ref. [22] and we compare 
the results with the description pointed in Ref. [19] based on the extension of (1.1). 

2.2 The stationary solutions and their stabihty 

The non zero stationary solutions of (j2.3p are solutions of the following second-order algebraic 
equation in the unknown L: 



PkL^ 



A 



-v'^ + Bin 



L + 



B2— + -4v^n 



0, 



(2.5) 



where we have put B2 = B1V2 and B^ = Aivi. 

Looking at the experimental results of Ref. [22], one notes that L is almost independent 
of V for V < Vc2, with a step change around Vd, while there is a variation of the slope near 
Vc2- We will concentrate on the change near Vc2- Reasoning as in Ref. [19], we observe that, 
under the hypothesis: 



^2 _ -B3 /5i ^ B3 

T~a;\J~a; 



(2.6) 



it follows that: 






'1 



Vl -V2 

and the solutions of equation (j2.5p can be written: 
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(2.7) 



(2. 



(2.9) 



El 

In the plane (F^, L), (j2.8p and (j2.9p represent two families of straight lines plotted in Fig. [H 
the first of them (equation (j2.8p l parallel to the V'^ axis and the second one (equation (|2.9p ) 
with a slope independent of fi. A linear stability analysis of these solutions shows that the 
solution (12.81) is stable if V is lower than: 



^c2 — ■ 



L 

Ai 



2— - — 

Ai /3. 



VLk, 



17 



Qk, 



(2.10) 



fl - V2 

(corresponding to the point of interception of the two straight lines (j2.8p and (j2.9p l. while, for 
values of Y higher than 1^2, the solution (|2.9p is stable. Therefore V"c2 represents the second 



critical counterflow-rotation velocity observed in the experiments of Ref. [22]. As we see, this 
critical velocity scales as ^/Ti, in agreement with experimental observations. 

The experimental data on the steady states of L allow us to determine the values assumed 
by the dimensionless quantities appearing in equation (12. Sh . One obtains: 



-4- = 0.0125, -4- = 3.90, -4- = 3.79, -^ = 0.025, (2.11) 

/3 (5 (5 (3 



from which we obtain: 



fi = 2.036, U2 = 0.97. (2.12) 



The coefficient /3, which controls the rate of evolution of L, cannot be determined from the 
knowledge of the stationary solutions. Comparing the value Ai/l3 obtained in the combined 
situation with the value (|2.4p obtained in absence of rotation, we can deduce that the first one 
is approximately 12 times the second one, which means that the coefficient Ai/P depends on 
the anisotropy and the polarity of the tangle, which depend on fi and V. 

Using the obtained values of the dimensionless quantities (|2.1ip and (|2.12p . the steady 
stationary solutions Lf and L2 become: 



Q T/2 Q 

Lf = 2.036-, and L^ = 0.0125^5- + 1.86- . (2.13) 

In Fig. [1] a comparison of such stationary solutions Lf and Lg with the experimental data of 
Swanson et al. is shown [22]. The conclusion of such a fit is that the stationary vortex line 
density L^ and L2, solutions of the alternative Vinen's equation in the combined situation, 
are in good agreement with experimental data of Swanson et al. 

In Ref. [19], the stationary solutions of equation (|2.2p had the form: 

i/2_,.ov./^ „.. ri/2 -^-^-y ^^- 1^ 



L'=l.i27\- and L'^ = 0.047— + 1.25W - (2.14) 

V K K \ K 

and the comparison with the experimental data led also to the conclusion that (j2.14p agree 
with the experiments by Swanson et al. [22J. Through (|2.13p and (|2.14p have a different 
mathematical form, in the range of the available experimental data, both of them lead to 
reasonable results. 

From such conclusions an interesting problem is to establish which equation, either (|2.2p 
based on the usual Vinen's equation or (j2.3p based on the alternative Vinen's equation, fits 
better the experimental data obtained by Swanson, Donnelly and Barenghi |22j . 

From a first comparison, the two stationary solutions (j2.14h ) and (j2.13b ) represent the 
same straight line in the plane (L, V'^) in the range V^^ < V^ < V^|. So, an eventual difference 
between both equations could be found in the range V'^ > V^. To do that, we calculate the 
errors a between L2 and the corresponding experimental value L, and a between L^ and L, 
respectively, in such a way that we can compare the accuracy of the two models. 

To find these errors, we consider the experimental values V^^ and J7j of the experiments to 
which L2ij, L2ij and Lij correspond, obtaining: 

a = J?l.li^^p^ = 963 (2.15) 



and 



Ei,j L^ij - h 



a^ = ^ =^±^ ^ = 419, (2.16) 

where N is the number of experimental data, which is equal for both cases. Prom (|2.15|) and 
(j2.16p we can establish that the stationary solution of the alternative Vinen's equation ap- 
proaches better the experimental data (for V'^ > V^) than that of the usual Vinen's equation. 

2.3 The first critical velocity 

The model based on the equation (|2.3p does not describe the existence of the first critical 
velocity Vd mentioned in Section (j2.ip . in which the value of L has a small steep change. To 
do this, we assume that the coefficient vi depends on fi and V as: 



1^1= All- B tanh 



N'\f,-C 



(2.17) 



with A, B and C constants. Thus, for V"^ <C V^i = -^kQ, it results v\ "^ A — BA and for 
V'^ » y^l, ui = A + BA, while the constant C is related to Vd, and 2B gives the size of the 
step of ui near Vd- In fact, if Vd is small, the domain of V in which the mentioned transition 
occurs is very narrow, as observed in experiments. 

In (|2.17p , the critical value Vd of the counterflow velocity is given in terms of coefficient 
Cby: 

K'l = ^^^- (2.18) 

Using the experimental values of Vd {Vd = 0.053\/l7 cm sec^^'^), it is seen that C = 
k/(0.053)2 = 0.355. 

To determine the coefficients A and B in (I2.17p . we consider that for a given value of 
il, for small values of V, the tangle will be completely oriented along the rotation axis, and 
vi = A — BA = 2. On the other side, when V 3> Vd (i-e. near Vc2) i^i assumes its higher 
value u'^"'^ furnished hy A + BA = y'^'^^ ; the value of u^^"'^ was obtained in Ref . [19] using 
experimental data of Ref. [22], and is i/™-"^ = 2L{Vc2)/Lr = 2.036. It is seen that the step 
in vi is indeed small. In (j2.17p . N' is a phenomenological coefficient characterizing the rate 
of growth of L near Vd and the experimental data show that N' > 20, but do not allow to 
determine it with precision. Here, we will chose for it the value N' = 22 proposed by Tsubota 
et al. in Ref. [18]. 

Expression (j2.17p is similar to that proposed in Ref. [19j and it is founded on the micro- 
scopic ideas about the nature of the transition, already proposed by Donnelly [5], according 
to which, for small V , the vortex lines are straight lines parallel to the rotation axis, but 
increasing values of V produce helical perturbations of the vortex lines around their low-^ 
configuration. The situation has been compared by Donnelly to magnetic systems, where the 
external field H contributes to the orientation of magnetic dipoles, while the temperature T 
has a disordering effect. Thus, the "tanh" term in ()2.17p is analogous to the expression de- 
scribing magnetization in terms of magnetic field and temperature in a 1/2 spin paramagnetic 
system. Other expressions, as for instance, Langevin's one for classical paramagnetism, could 
also be used [Ml |28l [29]. This ansatz is similar to that proposed in Ref. [30], to explain the 
transition from the laminar to the turbulent regime in pure counterflow, and is based on an 
analogous physical basis: there, the flow was producing the helicoidal excitation waves along 



the vortices pinned to the wahs of the container always present in the laminar regime. We 
recall that in pure counterflow, in containers with circular and square section, there are three 
different regimes distinguished by two critical counterflow velocities: under the first critical 
velocity we have the laminar regime where only a few of vortices pinned to the walls of the 
containers are present, between the two critical velocities a state of low vortex line density (TI 
regime) is formed, and, at last, above the second critical velocity a state of higher values of L 
(Til regime) is present. In Ref. [30] the transition from the laminar regime to the turbulent 
TI regime was explained supposing that small localized arrays of quantized vortices appear 
when the counterflow velocity reaches the first critical velocity, because Kelvin waves may 
be propagated in these pinned vortices. Whereas, when the counterflow reaches the second 
critical velocity, the TI turbulent regime, which is an inhomogeneous and locally polarized 
state, becomes unstable, with a transition to an homogeneous slightly not isotropic state TIL 
In the microscopic model we have commented on, the second critical velocity Vc2 is inter- 
preted as the velocity where the helical vortex lines produced in Vd have reached an amplitude 
of the order of the average vortex separation and have broken and reconnected, and form a 
disordered tangle. This explanation is analogous to the one given in Ref. [30] to explain the 
transition from TI to Til turbulent regimes. 



2.4 Non- stationary solutions of the Vinen's and alternative Vinen's equa- 
tions 



In this Subsection we study the non-stationary behavior of (12. 2p and (j2.3p . Though the lack of 
experiments about the evolution of the vortex line density L in this more general case (rotation 
and counterflow) does not allow us to compare directly our results with experimental data, 
however we can arrive at some interesting conclusions concerning the difference of behavior. 

First of all, we have to state that the analysis below refers to L as dependent variable 
when the growth of V or of is very small. Two main situations are considered, in the first 
one the angular velocity is fixed and the counterflow velocity moves between two consecutive 
experimental values (see Fig. [T|); in the second one the opposite situation is assumed, that is 

V is fixed and fi grows in a small range. This choice is due to the fact that when the ranges of 

V and r2 are sufficiently large, the coefficients of the Vinen's equations may be not constant, 
as showed by Schwarz and Rozen in Ref. [31], because they may depend on the anisotropy and 
polarization of the tangle. 

Denoting with Lq the initial value of L, the solution of the evolution equation ()2.3p is: 



l3Kt{L 
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(2.19) 



with Li and L2 given by (|2.13p : while the solution of the equation (j2.2p can be written as: 







(2.20) 



where Lq is the initial value of L, and Li and L2 are expressed by ()2.14p . 



In order to compare the unsteady solutions (|2.19p and (I2.20p of (12. 3p and (|2.2p . a value 
for the coefficient /? must be chosen. As already said, /? may depend on the anisotropy and 
polarization of the tangle, therefore it may have a different value with respect to the one in 
pure counterflow situation. However, since this dependence is not known in this section, to 
perform this comparison we choose the value of (3 in pure counterflow, namely (3 = 1/2tt. 

For Q. fixed. Now, we choose some values for V"^ and Q. in order to plot the solutions of the 
two models. First of all we consider the case V"^ < V'^21 ^'^d in particular the values V"^ = 0.0072 
and ri/27r = 0.4 to which the following values of the stationary solutions correspond: 

Li = L^ = 5132, L2 = 4455, L^ = 4779. 

For the initial value Lq we choose Lq = L/j = l^j k. Here, all the values for L, V"^ and Jl will 
be expressed in cm~^, cm^ s~^ and rad s~^, respectively. 

Prom the analysis of Ref. [19j we already know that in this range the stationary solution L\ 
is stable. The same conclusion is reached by looking at the plot of the non-stationary solutions 
(j2.19p and (j2.20p of the two models in Fig. [2j Further, we note that the values of L2 and Lg 
are smaller than L\ and that the non-stationary solutions approach to the stable stationary 
one, Li, in relatively similar times. 

Following the same process as above and setting the same value for $7 and a value V'^ = 
0.0626 slightly higher than V^2) we find the following values for the stationary solutions: 

Li = L^ = 5132, L2 = 5553, L^ = 5476, 
and for Ln two different values Lo, , and L9 are chosen respectively for the two 

^ ^|V-2=0.0482 ^11/2=0.0482 ^ -^ 

solutions (j2.19p and (j2.20p (see Fig. [3])). Note that in this case the value of Li is smaller 
than L2 and L2 • As we know from previous studies, in this range the stationary solutions L2 
and L2 for the Vinen's equation and alternative Vinen's one are stable. This is confirmed in 
Fig. H where ([239]) and (12:20]) are plotted. 

In Fig. [3] we also note a different behavior with respect to that in Fig. [2l in fact, the 
two non-stationary solutions L approach the corresponding stationary values Lg and L2 in 
rather different times with a ratio of about 1:3, respectively. So, the solution of the alternative 
Vinen's equation is faster than that of the Vinen's equation. 

Furthermore, if we plot the non-stationary solutions for a value of V'^ much higher than 
V^2! '^6 note that the ratio between the temporal scales is yet bigger than the factor 3. In 
fact, by setting the same value of i7 and taking V'^ = 0.1878, the corresponding values of the 
stationary solutions become: 

Li = L^ = 5132, L2 = 6910, L^ = 7050, 

and the graphics of the solutions (|2.19p and (|2.20p are shown in Fig. [H As initial data, we 
have chosen Ln = L2 o and Ln = L9 for the Vinen's equation and alternative 

Vinen's one, respectively. Looking at these unsteady solutions, we note that the solution of the 
alternative Vinen's equation approaches to Lg in a much shorter time than the other solution 
requires to approach L2, by a ratio of about 1:5. 

Note that the time scales in Fig. [3] [100-300 seconds] are much longer than those in Fig. [¥] 
[15-75 seconds]. This is not surprising because Fig. [3] corresponds to a situation which is much 
closer to the critical velocity Vc2 than that corresponding to Fig. U) Indeed, it is known that 
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the dynamics near critical points and phase transitions is much slower than in situations far 
from them. 

For V fixed. In the three situation considered before, the angular velocity is always 
constant {0,/2ir = 0.4) whereas the counterflow velocity increase from an initial value Vq to a 
final one V. The global behavior is the same when we suppose the opposite situation, that is the 
counterflow velocity is maintained constant and the angular velocity increases from an initial 
value Qq to a final one Q. In fact, by choosing the counterflow velocity V'^ = 0.024 < y^| and 
increasing O from Qq/2tt = 0.95 to S7/27r = 1 the following values for the stationary solutions 
are obtained: 

Li = L^ = 12831, La = 11345, L^ = 12023. 

For Lq the value of Li at 0,q/2tt = 0.95 is chosen. The plots of the unsteady solutions with the 
previous values is shown in Fig. [5l from which the same conclusion of Fig. [2] may be reached. 
In Fig. [6]the two non-stationary solutions (|2.19p and (|2.20p are shown when V'^ = 0.1626 > 
V^2 ^^'^ ^ increases from CIq/Itt = 0.95 to Q/ln = 1. There, by choosing the initial states 
Lq = L2|fj ,2^^(, gg and Lq = L^ _ , the two solutions L approach to the steady states L2 

and L2 : 

Li = L^ = 12831, L2 = 13967, L^ = 13767, 

in different times, in agreement with the previous situations. 

3 Conclusions 

The possibility of at least two reasonable evolution equations for the vortex line density L, 
namely (|l.ip and (|1.2p . was known since the early days in which Vinen proposed them. How- 
ever, detailed comparisons for them are very scarce [7l [9]. This was due, in part, to the fact 
that both of them lead to the same form for the steady state results, namely L ~ V'^, and that 
their unsteady solutions are not sufficiently different to reach a definitive conclusion on their 
relative merit. Here, we have carried out a detailed comparison of an extension of both equa- 
tions (jl.ip and ()1.2|) to the simultaneous presence of counterflow and rotation. The extension 
of (jl.ip was already studied in Ref. [19]. Here we have studied the analogous extension of 
(|1.2p . We have seen that in steady states the solutions of both equations, namely (|2.13p and 
(j2.14p have a different form but in the range of available experimental results both of them 
yield a satisfactory approximate description of the experimental data. However, a deeper com- 
parison of the experimental errors, in (j2.15p and (j2.16p . shows that the description based on 
the alternative Vinen's equation is slightly better than the one based on the most well-known 
Vinen's equation. 

A new aspect we have explored is the unsteady behavior of the solutions of these equations. 
Here, both equations exhibit remarkable differences, and we show that the solutions of the 
alternative Vinen's tend much faster to their steady-state values. In fact, this difference 
depends on the value of the counterflow velocity. For V^ = 0.0626, slightly higher than the 
critical velocity 1/|, the time required to reach the steady state solutions is 3 times shorter in 
the alternative Vinen's equation than in the usual Vinen's equation, whereas for V'^ = 0.1878 
the difference is still more remarkable, the time scale of the alternative Vinen's equation being 
5 times shorter than that for the usual one. Though we lack detailed experimental data on 
this unsteady behavior, we know that the time required to reach the steady state was less than 
10 minutes according to Swanson et al, when the counterflow velocity V is slightly above the 
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critical velocity Vc2 and it increases between two consecutive experimental values (see pag. 
191, Ref. [22] )• According to the results of the Fig. [3l the temporal scale of the solution of the 
usual Vinen's equation is closer to the observations than the temporal scale corresponding to 
the alternative equation, which tends too fast to the final result. Thus, it seems that the usual 
equation is preferable on these grounds. 

However, it must be stressed that the value of /3 used in our analysis has been /3 ~ 0.16, 
the value corresponding to pure counterflow, but the value of /3 could depend on the polarity 
of the tangle, as mentioned below expressions (I2.12p . The temporal scale of the solution of 
the generalized usual Vinen's equation could be set equal to the experimental value by setting 
P = 0.11; instead, to adjust the temporal behaviour of the generalized alternative Vinen's 
equation a much more radical change in the value of P should be made, setting /? = 0.03. 
However, such a drastic reduction of the value of P seems at odds with the fact that Ai/P in 
the presence of rotation is smaller than in pure counterflow, as mentioned below (j2.12p . Thus, 
the dynamical behavior of the usual Vinen's production terms seems more suitable than the 
modified one. Our work makes also evident the need of more detailed studies of the dependence 
of /? — and other coefficients — on the polarization of the tangle. 
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Figure 1: Comparison of the solutions (|2.13|) (continuous line) with the experimental data by 
Swanson et al. I22l. 



Figure 2: Evolution of the vortex line density L towards its steady state value for the general- 
izations of the usual Vinen's equation (j2.2p [dotted line] and the alternative Vinen's equation 
()2.3p [continuous line] for i7/27r = 0.4 and V'^ = 0.0072, lower than the critical value Vc2- 



14 



Figure 3: As in Fig. [21 but for r2/27r = 0.4 and V'^ = 0.0626 slightly above the critical value 
Vc2- Note that the steady solutions differ only in a 0.75%, whereas the difference in the time 
necessary to reach the steady state differs in more than 300%. The values of Lq are the 
unsteady solutions L2 and L2 at the same f2 and V^ = 0.0482. 
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Figure 4: As in Fig. O but for Q/2tt = 0.4 and V'^ = 0.1878 much higher than the critical 
value Vc2- The times necessary to reach the steady state differ in a 500% whereas the steady 
state values differ only in a 12.5%. The values of Lq are the unsteady solutions L2 and Lj at 
the same Cl and V"^ = 0.1626. 
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Figure 5: As in Fig. [2l for i7/27r = 1 and V"^ = 0.024 < 1/|, lower than the critical value Vc2- 
Here, the initial data Lq is the value of Li at il/ln = 0.95 and at the same V. 
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Figure 6: As in Fig. [2l for r2/27r = 1 and V'^ = 0.1626, above the the critical value Vc2- The 
values of Lq are the unsteady solutions L2 and L2 at the same V"^ and i7/27r = 0.95. 
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